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Abstract: The knee phenomenon of the cosmic ray spectrum, which plays an important role in studying the 
acceleration mechanism of cosmic rays, is still an unsolved mystery. We try to reconcile the knee spectra measured 
by ARGO-YBJ and Tibet-///. A simple broken power-law model fails to explain the experimental data. Therefore a 
modihed broken power-law model with non-linear acceleration effects is adopted, which can describe the sharp knee 
structure. This model predicts that heavy elements dominate at the knee. 
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1 Introduction 


The energy spectrum of cosmic rays, an important 
key to realize the origin and acceleration mechanism of 
cosmic rays, can be simply described by a power law over 
many magnitudes of energy froirr 10® eV to 10®° eV [ij, 
except for a few distinctive structures, such as the knee 
phenomenon around 4 x 10^® eV. The spectral power in¬ 
dex rapidly steepens from about —2.7 before the knee to 
—3.1 over the knee Q. Numerous works on cosmic ray 
propagation and acceleration mechanisms have 

been done to investigate the origin of the knee. These 
include the diffusive shock acceleration (DSA) spectra 
origin of the knee based on the theory of non-linear 
diffusive particle acceleration by shock waves from su¬ 
pernova remnants (SNRs) and the contribution 

of nearby pulsar wind producing very hard spectra in 
the power-law spectra in the knee region ^ 141. How¬ 
ever, the origin of the knee structure is still unconfirmed. 
Currently, the interpretations of this phenomenon are 
based on phenomenological models from experimental 
measurements. Many results have been given by direct 
observations with balloon detectors [Ml and indirect 
observations with air-shower detectors on high altitude 
mountains 3) 21, 22|. The CREAM experiment, as a 
balloon detector, has announced precise measurements of 
energy spectra for individual nuclei ranging from 2 TeV 
to 200 TeV [i^. As an air-shower ground-based detec¬ 
tor located 4300 m above sea level, the Tibet-/// array 
has presented an all-particle energy spectrum of primary 
cosmic rays from 100 TeV to 100 PeV, revealing a sharp 
knee structure around 4 PeV. Another air-shower exper¬ 
iment with the same altitude as Tibet-///, ARGO-YBJ, 


recently gave a new measurement of the energy spec¬ 
trum of hydro gen and helium nuclei, exhibiting a clear 
knee structure 241. 

The latest measurements released by ARGO-YBJ are 
investigated by the joint operation of Resistive Plate 
Ghambers (RPC) detectors and a Cherenkov telescope. 
This hybrid detection bridges the gap between balloon 
detectors and ground-based experiments, improves the 
shower energy resolution and enhances the capability to 
discriminate showers induced by light nuclei from events 
initiated by heavier nuclei. They yield clear evidence for 
a knee-like structure in the spectrum of hydrogen and 
helium nuclei under 1 PeV |24| . 

In the following discussions, we explain the knee spec¬ 
tra measured by ARGO-YBJ and Tibet-/// by phe¬ 
nomenological models. The superposition of energy spec¬ 
trum of individual components described by a simple 
broken power-law cannot reconcile the light components 
and all-particle energy spectra. We therefore propose a 
parametric model to reconcile them. In Sec. 2 we give 
two different scenarios and parameters to fit the experi¬ 
mental data. Sec. 3 presents the average mass. Finally, 
Sec.4 gives some discussion and conclusions. 


2 Energy spectra formulation 

Cosmic rays emitted from the sources are most likely 
accelerated in the strong shock fronts of SNRs by the 
DSA mechanism. The particles are deflected by chaotic 
magnetic fields, cross shock fronts frequently, and there¬ 
fore gain energy up to the PeV region. This acceleration 
leads to the observed approximate power-law spectrum. 
The energy spectrum is modified during diffusive prop- 
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agation. The spectra near the sources vary from that 
observed on the earth, which may due to nuclear spalla¬ 
tion or decay, ionization losses, leaka ge fr om the galaxy 
and solar modulation for low energies |25l427l | . We try to 
adopt two different phenomenological models to explain 
the observed energy spectra around knee structure. One 
is an expression of a simple broken power-law and the 
other is the broken power-law with non-linear modifica¬ 
tion term. 


2.1 Broken po-wer-law model and parameters 


The all-particle primary energy spectra derived from 
extensive air shower experiments can be described by a 
broken power-law function. We therefore use this kind of 
function to parameterize the observed differential energy 
spectra for individual elements of cosmic ray, which can 
be written as 


A. 

dE 


ZCb 


( 1 ) 


where ^ is the differential flux of the individual ele¬ 
ment with energy E, jo is the normalization constant, Cf, 
represents the break point of the proton spectrum, ^ is 
the atomic number, 7 is the power index in the energy 
range E « zeh, and 7 -I-A 7 denotes the power index 
when E » zCb- We try to use this expression to fit 
both the combined energy spectrum of hydrogen and he¬ 
lium nuclei observed by ARGO-YBJ and the all-particle 
energy spectrum given by Tibet-///. The interaction 
model SIBYLL for Tibet data is used in this work. We 
firstly fix the energy spectra of each component below the 
knee region, and then extrapolate the spectra beyond the 
knee. The observed data of individual elements (//, He, 
C, O, Ne, Mg, Si and Fe) from CREAM are adopted. 
For the spectrum index of each element, see [H, [^. An 
expression of simple power-law dj/dE = jE~^ is used 
to fit the spectrum of each component, with the fitting 
parameters j and spectral indices listed in Table 1. 

Table 1. Constant j and spectrum index 7 below 
the knee region. 


z 

Element 

j{eV ^sr 

7 

1 

H 

6.62x101® 

2.66 

2 

He 

6.40x101'' 

2.58 

6 

C 

2.43x101' 

2.61 

8 

0 

1 . 88 x 101 ® 

2.67 

10 

Ne 

1 . 68 x 101 ® 

2.72 

12 

Mg 

4.56x101' 

2.66 

14 

Si 

7.99x101' 

2.67 

26 

Fe 

6.71x101' 

2.63 


With Equation o, we fit the combined spectrum 
of hydrogen and helium nuclei observed by ARGO-YBJ 


18l . I24 l29j, and get the fitting result of energy break 
point £{, = 4 X and A 7 = 0.68 for the proton. 

The break point for other elements can be calculated by 
ztb- The spectral indices of all elements at the higher 
energy range beyond the knee region are assumed to 
be the same value as the result from ARGO-YBJ, i.e. 
—7 — A 7 = —3.34. We sum up the energy spectra of 
all elements to get the approximate all-particle spec¬ 
trum, seen as the red solid line in Fig. [T] It can be 
seen from Fig. [l]that this all-particle spectrum based on 
ARGO-YBJ observation with this simple broken power- 
law model cannot match the observed all-particle spec¬ 
trum given by Tibet-///. The simple broken power-law 
function therefore cannot describe the energy spectrum 
in the knee region, which perhaps is due to the influence 
of the DSA mechanism. 



E(eV) 


Fig. 1. All-particle spectrum (red solid line) cal¬ 
culated as the sum of individual components, 
with spectra of some individual elements. The 
data are cited from CREAMjp,!!^, ARGO-YBJ 

El, [13, [H and Tibet-/// [ilL 


2.2 Non-linear model and parameters 


Since the differential energy spectrum in the form of 
the simple broken power-law cannot reconcile the obser¬ 
vations of ARGO-YBJ and Tibet-///, a new correction 
term is added, so the broken power-law can now be ex¬ 
pressed as 


dE 


= joE--[l + a{—f](l + —)-^\ 
ZCb ZCb 


( 2 ) 


where a is assumed to be a z-dependent parameter, and 
,5 is a constant. The DSA mechanism may result in a 
spectrum with concave-up curvature, so we multiply the 
non-linear term l-|-a(^)^ into the simple broken power- 
law to structure a concave-up form as a phenomenolog¬ 
ical model. Without this non-linear term, the superim¬ 
posed all-particle spectrum from all elements is smooth 
and cannot fit the sharp knee structure as measured by 
Tibet-///. In addition, we have tried some parameters 
in a non-z-dependent form, but these cannot match the 
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data, so we finally adopt the form in Equation ([2|). An 
equation to express the non-linear effects is proposed in 
[SOj. They propose the non-linear process by a modi¬ 
fied exponential cut-off expression for the spectra from 
sources, and then deduce the observed spectra by the 
superposition of these spectra. From their formula, the 
model they deduced is a kind of broken power-law multi¬ 
plied by a complex modified term. While our expression 
form is simpler, we modify the simple broken power-law 
directly as the observed spectra and then parameterize 
this scenario. Combined with the newly published mea¬ 
surements from ARGO-YBJ, our model is easier to quan¬ 
tify the parameters precisely. A model with precisely 
confirmed parameters will help us to review its intrinsic 
physical issues. From Equation ([2|), we get a good fit 
result to match the observations from both ARGO-YBJ 
and Tibet-/// with parameters a = 1.50-1-0.352, /3 = 1.20, 
£(, = 2.0 X lO'^^eV and A 7 = 1.88 for the proton. Ay val¬ 
ues of other elements are assumed to obey the formula 
—7 —A 7 -I -/3 = —3.34. The fitting spectra are plotted in 
Fig. [21 where the solid red line represents the all-particle 
spectrum. From Fig. [2] we can see that the all-particle 
spectrum based on ARGO-YBJ observation can match 
the observed all-particle spectrum given by Tibet-///. 
Based on this model, the distinctive knee structure is 
mostly contributed by heavy elements such as Fe. 



E(eV) 


Fig. 2. All-particle spectrum (solid red line) cal¬ 
culated as the sum of individual components, 
with spectra of some individual elements. The 
data are cited from CREA M [^ . [2^ . ARGO-YBJ 
[H, m, [ 2 ^ and Tibet-/// l2l|?^ 


In order to obtain the best fitting parameters for the 
non-linear model, we employ the Markov Chain Monte 
Carlo (MCMC) [1^ 01: a very efficient algorithm, to 
achieve our goal. Considering that the calculation of 
will be affected by the too small errors of all-particle 
spectrum given by Tibet-///, we amplify the errors of 
this spectrum by taking account of energy resolution. 
The initial form of data with errors listed in the liter¬ 
ature 2l| are expressed by /" ± AF, where F denotes 


the flux and AF denotes the systematic error. As a new 
deduced additional error, a new term is added to this 
expression, i.e. F±AF±F{'y — l){AE/E), where 7 de¬ 
notes the spectrum index. We take the energy resolution 
AE/E = 0.15 as a constant. The following 5 parameters: 
a =p]_+p 2 Z, /3, £6 and Ay, are set free for the MCMC 
algorithm, where a consists of free parameters pi , p 2 and 
redshift z. The parameters £5 and Ay here refer to pro¬ 
tons as mentioned above. Because every spectrum index 
of elements in the higher energy band beyond the knee 
is supposed to be the same value, we can get the Ay of 
other components from the Ay of the proton. By this 
MCMC algorithm, we obtain the best fitting parameters 
as follows: a = 2.21 -f 0.452, P = 1.62, £5 = 1.0 x W^eV 
and Ay = 2.25 for the proton. The parameters a, (3 and 
£{, decide the shapes of the energy spectra. Different val¬ 
ues result in different spectra, as shown in Fig. |3l where 
the elements // and Ee are taken as examples. 






Fig. 3. Energy spectrum shapes determined by a, 
/3 and tb for elements H (left) and Fe (right). The 
upper two figures represent the dependence of the 
spectra on a, the middle two figures represent the 
dependence on /3, and the lower two figures repre¬ 
sent the dependence on £ 5 . The solid red lines in 
all six figures represent the currently used values. 


The best fit result is shown in Fig. |4l The measure¬ 
ments of hydrogen and helium nuclei spectra using the 
SIBYLL model by the KASCADE group [sij are com¬ 
pared with our model in Fig. U) The KASKADE group 
claimed that the energy spectra of both proton and he¬ 
lium show a knee-like feature, but in our model, this 
seems not to be the case. 
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E(eV) 


Fig. 4. All-particle spectrum (solid red line) cal¬ 
culated as the sum of individual components, 
with spectra of some individual elements. The 
data are cited from CREAM ARGO-YBJ 

[ii,[2i,0, Tibet-/// 0 and KASCADE 


3 Average mass 

A commonly-used quantity to characterize the mass 
composition of cosmic rays is the mean logarithmic mass, 
which is defined as 

{InA) = '^rilnAi, (3) 


Ti being the relative fraction of nuclei of mass Ai. The av¬ 
erage mass is reported by many air shower experiments, 
but the results are divergent because of the poor primary 
mass resolution. We calculate the average mass for two 
cases: one for the former non-linear model and another 
for the latter best fitting case, as shown in Fig. [S] 

From our derived average mass of cosmic rays, we can 
see {InA) goes up between several hundreds of TeV and 
several tens of PeV. It implies that heavier components 
occupy an increasing proportion with the increase of en¬ 
ergy, and predicts an iron-dominant composition above 
the knee. But the components above the knee are not 
all Fe , because {InA) for Fe is about 4.0. Further ex¬ 
perimental measurements of the chemical composition 
of high energy cosmic rays will verify our prediction. If 
the light components dominate the knee structure, this 
model can be excluded, while if the heavy components 
dominate the knee structure, this model will be con¬ 
firmed, and the measurements will be helpful to improve 
the model parameters. 



Fig. 5. Average mass calculated for the non-linear 
model, compared with that of Tibet-/// (blue 
dashed line). The green dashed line represents 
the first fitting result by the non-linear model 
and the red solid line represents the best fitting 
result. The data are cited from RUNJOB [3^ . 
JACEE [13, CASA/MIA [H^AS A/BLANC A 
[li], KASCADE 27thICRC H, BASJE-MAS 
and BASJE-Cherenkov [j^l- and Fly’s Eye 

4 Discussion and conclusion 

Many mechanisms have been discussed to explain the 
knee structure of cosmic rays, including the diffusive 
shock acceleration mechanism of cosmic rays in super¬ 
nova explosion shock fronts, leakage from the galaxy, in¬ 
teractions with background particles in the interstellar 
medium, and new physics of interactions in the atmo¬ 
sphere. Some models are likely to be excluded, however, 
such as the leaky box model, interactions with back¬ 
ground particles and new types of interactions in the 
atmosphere. The leakage from the galaxy does not give 
a distinct change of the spectral slope at the knee. The 
scenario of interactions with background particles yields 
a very light mean logarithmic mass, while the measure¬ 
ments indicate an increase of mass with energy. For 
new types of interactions in the atmosphere, there are no 
inconsistencies between the different air shower compo¬ 
nents by the simultaneous observations from KASCADE, 
and the measurements of different air shower components 
can be interpreted with standard particle physics (38| - 
We propose two scenarios to reconcile the light com¬ 
ponents and all-particle energy spectra by phenomeno¬ 
logical models. In the first scenario, we use a simple 
broken power-law model to fit the combined energy spec¬ 
trum of hydrogen and helium nuclei measured by ARGO- 
YBJ and the all-particle spectrum measured by Tibet- 
///, but the all-particle spectrum cannot be fitted by 
this model. Hence, the simple broken power-law fails to 
explain the knee structure. In the second scenario, a non¬ 
linear modification term is added to the broken power- 
law. This modified broken power-law model with non- 
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linear acceleration effects can describe the sharp knee 
structure. The mechanism of non-linear effects at shock 
fronts of SNRs has been explained in earlier works [s^. 
Particles with higher energies can scatter farther ahead 
of the shock, which results in an effectively higher com¬ 
pression ratio, so a more locally hard spectrum is pro¬ 
duced, and a concave-up curvature to the accelerated- 
particle distribution is expected (40l |. In addition, this 
modified model predicts that heavy elements dominate 
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